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|> ; ABSTRACT 

Phase noise caused by an inhomogeneous, time- variable water vapor distribu- 
tion in our atmosphere reduces the angular resolution, visibility amplitude and 
coherence time of millimeter and submillimeter wavelength interferometers. We 
present early results from our total power radiometry phase correction experiment 
carried out with the Submillimeter Array on Mauna Kea. From accurate mea- 
surements of the atmospheric emission along the lines of sight of two elements 
Q-V of the array, we estimated the differential atmospheric electrical path between 

them. In one test, presented here, the phase correction technique reduced the 
rms phase noise at 230 GHz from 72° to 27° over a 20 minute period with a 
2.5 second integration time. This corresponds to a residual differential electrical 
path of 98 /xm, or 15 /zm of precipitable water vapor, and raises the coherence in 
^ ' the 20 minute period from 0.45 to 0.9. 

. ?? : 

Subject headings: atmospheric effects — instrumentation: adaptive optics — site 
testing — techniques: interferometric — submillimeter 



1. Introduction 

The water vapor in our atmosphere is spatially inhomogeneous and highly time- variable, 
resulting in atmospheric electrical path fluctuations that alter the phase of propagating elec- 
tromagnetic waves. At submillimeter wavelengths, the effect of atmospheric fluctuations is 
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particularly severe, and limits the resolution and coherence of interferometric arrays. With- 
out atmospheric phase correction schemes, submillimeter wave interferometers such as the 
Submillimeter Array 1 (SMA) and the upcoming Atacama Large Millimeter Array (ALMA) 
will not be able to image at their diffraction limits in cases where self-calibration or fast 
switching is not possible. The advantages and limitations of fast switching are described by 
both Carilli, & Holdaway (1999) and Lay (1997b). At the SMA, fast switching may have 
limited value because of the large angular distances between calibrator sources of sufficient 
strength (Battat 2004). 

The atmospheric phase at millimeter and submillimeter wavelengths can be determined 
empirically by monitoring the continuum emission at the same frequency as the astronomical 
observation. The use of a single receiver for both source observation and phase calibration 
makes it possible to monitor the volume of atmosphere responsible for the phase errors. 
The accuracy of the phase correction is then independent of baseline length. Such a phase 
correction scheme can not only improve the quality of observations, but also broaden the 
range of weather conditions during which observations are possible. This translates to a 
more efficient observing schedule. 

We have developed highly stable total power detectors and have explored how reliably 
they can be used for accurate phase correction. The initial results described here show that 
total power radiometry has great promise to solve the problem of image blurring caused by 
atmospherically induced phase errors. 

2. Atmospheric Phase Errors 

Atmospherically induced phase errors degrade the measured visibility amplitudes, image 
resolution and dynamic range of submillimeter interferometers (e.g. Carilli, & Holdaway 
1999). Although water vapor is essentially non- dispersive at radio frequencies (one mm of 
precipitable water vapor, pwv, creates ~6.5 mm of excess electrical path), the corresponding 
phase errors increase linearly with frequency. In the immediate vicinity of strong atmospheric 
lines, the water vapor is dispersive, but at these frequencies the high opacity makes ground 
based radio astronomy impossible (Sutton, & Hueckstaedt 1996). In the well established 
Kolmogorov description of atmospheric turbulence, the root phase structure function, which 
describes how phase errors depend on the antenna separation or baseline, b, is a broken power 



1 The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and 
the Academia Sinica. For more on the SMA and its characteristics see Ho et al. (2004). 
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law proportional to 6 /3 / 2 , where (3/2 is called the Kolmogorov exponent. For baselines shorter 
than the thickness of the turbulent layer (b ^ 0.5 — 2 km), the turbulence is approximately 
isotropic in three dimensions and (3/2 = 5/6. For baselines longer than the turbulent layer 
thickness, the turbulence is two dimensional and (3/2 = 1/3 (Carilli, & Holdaway 1999). 
On size scales larger than the outer scale of atmospheric water vapor (typically 5-10 km), 
the phase errors are independent of baseline length. Because of this outer scale and the 
associated long time scales, Very Long Baseline Interferometry (b > 1000 km) is viable. 

The expected rms phase error between antennas pointed to the zenith can be written 

as 



where v is the observing frequency, 6 EL is the elevation angle of the observation and 7 is 
1/2 or 1 for three or two dimensional turbulence, respectively. Measurements at the SMA 
site, in the three dimensional turbulence regime, give K ~ 30° and (3/2 ~ 0.75, though 
large variations in both parameters are seen (Masson 1994). If the distribution of phase 
fluctuations is Gaussian, then the expectation of the measured visibility amplitude, {V m ) 
(where (. . .) denotes an ensemble average) is related to the true visibility amplitude, V, by 



(Thompson et al. 2001). Thus an rms phase error of 1 radian, which is expected on a 140 m 
baseline zenith observation at 345 GHz (Equation 1), introduces a 40% coherence loss. The 
resolution of the astronomical image is degraded by a convolution of the true image with a 
seeing function that depends on the root phase structure function (Thompson et al. 2001). 
This effect is the sub millimeter analog to optical seeing. Masson (1994) showed that the 
median resolution allowed by the Mauna Kea atmosphere is 1'.'2 at 345 GHz, corresponding 
to an effective maximum baseline of about 100 m. 

The dominant contribution to fluctuations in atmospheric delay comes from poorly 
mixed, wind-blown water vapor cells of various sizes that constitute a turbulent region con- 
fined to the first few kilometers above ground level (Woody et al. 2000). The relevant time 
scale for these fluctuations is ~ b/v where v is the wind speed along the baseline. At the 
SMA, b is 10 to 500 m and v is typically 5 to 10 m s _1 , so that the dominant fluctuations 
have periods of 1 to 100 seconds. 

Because surface meteorological measurements cannot provide an accurate determination 
of total water vapor content (Waters 1976), remote sensing measurements, such as microwave 
radiometry, are necessary. The real and imaginary parts of the dielectric constant of water 
vapor are related by the Kramers- Kronig relation and thus there is a strong correlation 




(1) 



{Vm) 




(2) 



V 



-4- 



between atmospheric opacity and phase delay (Thompson et al. 2001). The opacity in the 
atmospheric continuum windows below 1 THz is dominated by the far wings of infrared 
water resonance lines and also includes a continuum component that is proportional to the 
square of frequency. The correlation between the atmospheric brightness temperature and 
water vapor content has been established empirically by comparing radiometer measurements 
with radiosonde data (e.g. Staelin 1966; Westwater, & Guiraud 1980; Moran, & Rosen 1981; 
Guiraud et al. 1979). 

3. Past Results 

The correlation between atmospheric phase delay and atmospheric brightness tempera- 
ture has been measured at frequencies near water vapor resonance lines and in the continuum 
between the lines (e.g. Welch 1999). 

The Owens Valley Radio Observatory (OVRO), using 22 GHz water line radiometers, 
has demonstrated a reduction in the rms phase error over 25 minutes at 100 GHz from 64 
degrees to 19 degrees, corresponding to a residual excess electrical path of 160 /im (Woody 
et al. 2000), and have shown that phase correction improved their astronomical image (Mar- 
vel, & Woody 1998). Wiedner et al. (2001) used a 183 GHz water line monitor to correct 
for atmospheric phase errors in an observation at 354 GHz with the Caltech Submillimeter 
Observatory (CSO) and James Clerk Maxwell Telescope (JCMT) as a two element inter- 
ferometer of baseline 164 m. The rms phase over 30 minutes was improved from 60 to 26 
degrees (61 /im of residual excess path). 

Total power phase correction at the frequency of astronomical observation, which re- 
quires exceptional receiver gain stability of a few parts in 10,000, has also achieved some 
success. Zivanovic et al. (1995) used this technique at the Berkeley-Illinois-Maryland- 
Association Millimeter Array (BIMA). The total power difference at 88 GHz correlated 
with slow drifts in the measured phase but there was significant disagreement on short (~30 
second) and long (~30 minute) time scales. Bremer et al. (1996) at the Institut de Radio 
Astronomie Millimetrique Plateau de Bure Interferometer (IRAM, PdBI) reported a consis- 
tent correlation between the total power difference, measured at 90 and 230 GHz, and the 
interferometric phase at those frequencies, under clear sky conditions. They improved the 
rms of phase data from > 100 degrees to 25-35 degrees at 230 GHz over a 1 minute time 
scale (100 /im of path). In the presence of clouds, the total power correction technique is 
seriously compromised because of liquid and frozen water (Waters 1976). 
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4. Total Power Radiometry for Atmospheric 
Phase Correction at the SMA 

The SMA is an interferometric array on Mauna Kea with baselines up to 500 meters 
and frequency coverage from 180 GHz to 900 GHz. At the time of this experiment, integra- 
tion times, t, as small as 2.5 seconds were possible. Gain variations, caused by cryocooler 
temperature fluctuations, typically limit the Superconductor Insulator Superconductor (SIS) 
heterodyne receiver (Blundell et al. 1995) stability to 0.1-1%. 

On Mauna Kea, the median rms differential path on a 100 m baseline is about 112 /im 
(Masson 1994), which corresponds to atmospheric brightness temperature fluctuations of 
0.25 K at 230 GHz. Typical system temperatures at 230 GHz at the SMA are 125 K, so the 
atmospheric temperature fluctuations are 500 times less than the system noise temperature. 
Thus a receiver stability of order 1 part in 5,000 is desired to ensure accurate measurement 
of atmospheric temperature fluctuations. 

A servo system that monitors the physical temperature of the mixer and adjusts the 
receiver gain to stabilize the receiver output has been developed and tested (Battat et al. 
2004). Our two matched total power detectors operate near room temperature and are 
temperature stabilized to 1 part in 10 5 . The rms stability of the continuum detector, with 
a 33 ms integration time, was measured at 1 part in 8,000 over 10 minutes, which is only 
15% greater than the fundamental 1/VWt thermal noise. The servo system was inactive 
during the tests described in this paper. Instead we relied on the intrinsic stability of the 
two cryocoolers which are known to be very stable. The stability of the entire receiving 
system was measured at 5 to 7 parts in 10 4 on a fixed temperature load. 

Systematic instrumental noise due to radiometer gain fluctuations translates to a noise 
limit in the measurement of the differential electrical path between antennas, AL, that 
depends on the system noise temperature, T sys , and atmospheric conditions in the following 

way 

at ~o*f dL /dpwv \ ( dT bri /dpwv y l ( T sys \ ( AG/G \ 
AL ~ 35 { 6 .5 J^ieKmm-^ ^ 1251^500=* J ^ (3) 
where AG/G is the fractional receiver stability and T\, r i is the atmospheric brightness tem- 
perature which is approximately ToTocsc(# EL ) where To is the zenith opacity and To is the 
physical temperature of the atmosphere. Also, T sys « T rx + T^j where T rx is the receiver 
noise temperature. The value of dT bri j 'dpwv is computed from a mean atmospheric model for 
Mauna Kea in January (Paine 2004). At 230 GHz it is largely independent of water vapor 
content, ranging from 16.4 to 15.8 K mm" 1 for 1 to 4 mm pwv. Thus for AG/G = 1/2000 
in median 230 GHz observing conditions (dL/dpwv ~ 6.5, dT^ij 'dpwv ~ 16 K nun . T S y S ~ 
125 K) the systematic instrumental noise is 35 /xm, or 10 degrees of phase, for zenith obser- 
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vations. The thermal noise in the detector contributes only 1 /xm path error in 2.5 seconds 
of integration. 

5. Observations and Results 

We made test astronomical observations during several nights in January, 2004 with 
various integration times, observing frequencies and sources. Here we present 20 minutes of 
data on the quasar 3C273 at 230.5 GHz with two antennas on a 140 m north-south baseline 
from 16:35 UTC to 16:55 UTC with a 2.5 second integration time. Over this time the source 
azimuth and elevation spanned 246.7-250.7 and 51.7-47.0 degrees, respectively. The surface 
air temperature ranged from -5.1 to -4.7 °C, and the wind speed ranged from 5 to 15 m 
s _1 coming from the east and east-south-east. The 225 GHz zenith opacity was between 
0.13 and 0.14 during the first part of the observation and then rapidly increased to 0.18 
15 minutes into the observation. This opacity change was accompanied by the onset of 
large amplitude phase fluctuations (see Figure 1). Note that increased turbulence is not 
necessarily associated with increased total opacity (see e.g. Hinder 1972). The opacities 
listed above correspond to zenith pwv content between 3 and 4.2 mm, which is considered 
poor weather (75th percentile) on Mauna Kea (Masson 1994). At 1 part in 25, the rms to 
mean fluctuations in the individual total power signals were easily detectable above our 5 
parts in 10 4 noise floor. 

The top panel in Figure 1 shows S ma, the interferometer phase as measured by the 
SMA correlator and TP the total power phase prediction which was computed from the 
total power measurements made at two antennas, Pi and P 2 , by the equation 

TP = (AP 1 - BP 2 ) + C. (4) 

The scaling factors, A and B, were required because we did not have accurate temperature 
scales for either Pi or P 2 , and C was required because of the 2n ambiguity of the interfer- 
ometer phase, 4>sma- 

To compare </>sma with TP , we binned the total power data, sampled 
at 30 Hz, onto the 2.5 second SMA integrations. We determined A, B and C by minimizing 
the mean squared error of S ma — 0tp- We did not remove any temporal drifts from Pi, P 2 

Or S MA- 

Subtracting the total power phase model from SM a improves the rms error from 72° 
to 27°, corresponding to a residual differential electrical path of 98 /im or a differential 
water column of 15 /im. Figure 2 shows the strong correlation between TP and S ma- The 
coherence factor (see Equation 2) is improved from 0.45 to 0.90 by the total power phase 
correction. During the observation, the source elevation was ~50°, thus the noise floor was 
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11° (see Equation 3). The expected thermal phase noise of the SMA for 3C273 (with a flux 
of 4 Jy at A = 0.9 mm) was about 5°. These independent non-atmospheric phase errors 
add in quadrature to contribute ~12° to the residual phase. The remaining error in phase 
is attributed, in part, to unmodeled time variations in the A, B and C parameters and the 
lack of perfect correlation between total power and phase. 

The root temporal power spectrum, <j)(f), describes the rms phase error as a function 
of frequency, /. Kolmogorov theory, under the assumption that the turbulent velocity is 
small compared with the wind speed (frozen turbulence, or the Taylor hypothesis), predicts 
4>{f) oc f~ a with a = {3/2 + 0.5. Thus a = 4/3 above the corner frequency, f c ~ v/b, and 
a — 1 = 1/3 below f c (Masson 1994). The root temporal power spectrum of our observation is 
shown in Figure 3. The best fit function has a = 1.0 (slightly shallower than the Kolmogorov 
prediction) with a corner frequency of ~ 0.02 Hz. This corresponds to a time scale of 50 
seconds or a wind speed along the baseline of ~ 3 m s -1 , consistent with surface meteorology. 
This type of spectral break has been reported by others (e.g. Lay 1997a; Masson 1994). 

The root temporal power spectrum of the residual phase has a = 0.41, significantly 
flatter than the raw phase spectrum. With data from a single baseline, it is difficult to 
determine the residual root phase structure function from the residual temporal power spec- 
trum in order to extrapolate our results to longer baselines. However, a worst case scenario 
would be to assume that the phase errors continued to grow with baseline as (3/2 = 0.75 as 
is typical for uncorrected measurements on Mauna Kea. Then diffraction limited imaging 
would be possible on baselines as large as 380 m. It is more reasonable that {3/2 < 0.6, in 
which case diffraction limited imaging would be possible on all SMA baselines. 

6. Discussion and Conclusions 

The data presented above demonstrate that total power radiometry at the SMA can 
provide an accurate measure of the atmospheric phase between antennas and can help recover 
diffraction-limited image resolution. The structure seen in the residual phase (27° rms) is 
partly caused by instrumental noise from the total power radiometers and interferometer 
(see Section 5) and could also be due to changes in the vertical profile of the water vapor, 
or the presence of liquid water which contributes significantly to the atmospheric brightness 
temperature but not to the electrical path (Carilli, & Holdaway 1999). In fact, when the data 
set was divided into two 10 minute segments and A, B, C were derived for each segment, the 
resulting rms of the entire observation was improved to 22° rms, or 83 /xm of electrical path, 
which suggests that these parameters are sensitive to changing meteorological conditions. 
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In the future, the total power phase correction could be applied in real time by rapidly 
adjusting the phase of the local oscillator in each antenna, or it could be done during data 
analysis. In any case, the A, B and C coefficients in Equation 4 should be empirically 
determined from measurements of a phase calibration source and then applied to the science 
source observations that immediately follow. To achieve this goal, further studies must 
be undertaken to determine how rapidly A and B vary in time, and over what range of 
elevations they are valid before changes in contaminating factors such as ground pickup 
become significant. In addition, when multiple baselines are present in an observation, it 
will be possible to derive the A and B parameters from a simultaneous fit to the data from 
all baselines. This approach would provide a more robust parameter determination. 

The development of a consistently accurate system to remove atmospheric phase delay 
is of immediate importance to the development of submillimeter wave interferometry. At 
present no such system exists, but our data show that total power radiometry may be a 
viable option and is worthy of further development. 

We thank R. Christensen, T. R. Hunter, R. Kimberk, P. S. Leiker and M. Wiedner for 
their help with instrumentation. J.B. was supported by the NDSEG graduate fellowship. 
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Fig. 1. — Phase correction data from 2004 Jan 30 tracking the quasar 3C273 for 20 minutes 
at 230.5 GHz. Top: The solid line shows S ma, the phase as measured by the SMA with a 
2.5 second integration time, while the broken line shows TP , the phase determined from the 
total power measurements at the same frequency as the science observations. The opacity 
increased sharply after the 15 minute mark. Bottom: The residual phase, S ma— 0tp- 
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Fig. 2. — Correlation between the interferometer phase, S ma, (abscissa) and the total power 
phase, TP , as given by Equation 4 (ordinate). The solid line has unity slope and zero 
intercept. 
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Fig. 3. — Top: Root phase power spectrum of S ma (thick line) and S ma — 0tp (thin line). 
Also shown are the Kolmogorov prediction (a = 4/3, dotted line) and the best fits to the 
raw and residual data between frequencies marked with arrows (a = 1.0 and a = 0.41, 
respectively). Bottom: The ratio of the residual phase and S ma plotted on log-linear axes. 



